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■ Abstract: The theory of particle acceleration at shock fronts is briefly reviewed, 
with special emphasis on the production of the particles responsible for the non- 
thermal emission from blazars. The flat radio/IR spectra of these sources cannot 
be produced by diffusive acceleration at a simple nonrelativistic shock front prop- 

^ I agating in a homogeneous medium. It can, however, be produced by a single 

' unmodified mildly relativistic shock, if the pressure in the shocked gas is provided 

O . by the leptonic component, or, independently of the equation of state, by a rela- 

I tivistic shock which is oblique to the magnetic field. The analytic theory of these 

shocks makes several simplifications, but Monte-Carlo simulations exist which 
. extend the range of validity. Of particular interest is acceleration in a tangled 

magnetic field. Here, however, the Monte-Carlo simulations have not yet yielded 
Q^' unambiguous results. The 'homogeneous' models of blazar emission are discussed, 

Q , and is it shown that they imply a geometry of the emitting region which is laminar 

I in form, with an aspect ratio of d/R < 3% in the case of Mkn 421. Identifying 

■ these with relativistic shock fronts, a model of acceleration is described, which 
. . ', displays characteristic variations in the synchrotron spectral index with intensity. 
> ■ 

X 

1 Introduction 



The problem of how particles are accelerated to nonthermal energies in 
relativistic jets has been discussed for a number of years (e.g., Begelman, 
Blandford & Rees 1984). Diffusive acceleration at nonrelativistic shocks is 
a possibility which appears to provide a reasonable picture of acceleration 
at several jet hot-spots which emit optical synchrotron radiation (Meisen- 
heimer et al. 1989). However, the rapid variability of emission from blazar 
jets implies substantial doppler boosting, which suggests that the nonrel- 
ativistic theory may be inappropriate. Relativistic shocks have also been 
known for some time to be potentially effective accelerators (Kirk & Schnei- 
der 1987), but, unlike their nonrelativistic counterparts, they do not lead to 
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a unique spectral index which is independent of the details of the scatter- 
ing process and the shock speed. Instead, a wide range of spectral slopes is 
possible. 

The best developed model of blazar spectra is based on the picture of a 
shock front moving down a jet (Marscher & Gear 1985). Both the magnetic 
field and the particle distribution vary with position, and the observed radi- 
ation is a superposition of emission from different parts of the jet (Ghisellini, 
Maraschi & Treves 1985, BaUard et al. 1990, Hughes, AUer & AUer 1991, 
Maraschi, Ghisellini & Celotti 1992, Levinson Sz Blandford 1995, Marscher 
Kl Travis 1996). These models are generally called 'inhomogeneous'. How- 
ever, recent results concerning the rapid variability of blazar sources at all 
frequencies (Wagner & Witzcl 1995) and the detection of their emission 
at energies of up to at least 10 TeV (Aharonian et al. 1997) have provided 
new restrictions on the possible acceleration mechanisms. Especially the ob- 
served simultaneous variations in X-rays and TeV gamma-rays indicate that 
a single population of particles in a relatively localised region is responsible. 
Consequently, 'homogeneous' models have been widely discussed (Dermer & 
Schhckeiser 1993, Macomb et al. 1995, 1996, Bloom & Marscher 1996, Ghis- 
ellini, Maraschi & Dondi 1996, Inoue & Takahara 1996, Stecker, de Jager & 
Salamon 1996, Comastri et al. 1997, Mastichiadis & Kirk 1997) and there 
is rough agreement on the parameters of the emission region. 

Simultaneous variability of emission at widely differing wavelengths is 
more easily interpreted in terms of directly accelerated electrons than in 
the hadronic models (Mannheim 1993, Mannheim et al. 1996), in which the 
gamma-ray emission stems originally from ultra-relativistic protons, and 
only the optical and (in some sources) X-ray photons arise from directly 
accelerated electrons. Protons and electrons are accelerated in regions of 
different spatial extent on very different timescales, because the waves which 
provide the scattering centres have a wavelength comparable to the particle's 
gyro radius, and the acceleration time of a particle is probably somewhat 
longer than its gyro period. Consequently, one would expect the radiation 
produced by relatively low energy (GeV to TeV) electrons to have different 
variability properties to that produced by protons of energy above 1 EeV. 
However, further information concerning the emission of these objects at 
energies above 20 TeV (Meyer & Westerhoff 1996) may well provide the 
most reliable way distinguishing between the leptonic and hadronic models. 
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2 Basic properties of shock acceleration 

It is well-known that particles which diffuse in the neighbourhood of a non- 
relativistic shock front can be accelerated into a power-law distribution 
such that the density of particles with Lorentz factor 7 is a power-law, 

77.(7) 0^ 7"'^"''^; with 

S = 3pe/(Pc - 1) (1) 

where pc is the compression ratio of the shock front. The physical basis 
of acceleration is that an energetic particle scatters elastically off magnetic 
fluctuations in the background plasma, which enables it to cross and recross 
the shock front. Simple kinematics lead to an energy gain on each crossing, 
since, at a shock front, the scattering centres take part in the plasma com- 
pression. This competes with the possibility that a particle moves off into 
the downstream plasma and does not return, to yield a power-law spectrum 
(for a review see Kirk, Melrose & Priest 1994). 

Before they have had time to cool, the synchrotron emission of electrons 
accelerated into a distribution described by Eq. (1) is a power law I{u) oc 
with a = {s — 3)/2. Thus, a strong nonrelativistic shock front in a 
fully ionised gas, which has pc — 4, produces a = 0.5. There are several 
effects which modify this result. However, most of them lead to a softening 
of the predicted spectrum, at least for lower energy electrons (e.g. Bell 1987, 
Ellison & Reynolds 1991, Duffy, BaU & Kirk 1995). In the case of blazars, the 
homogeneous model requires a harder spectrum - for Mrk 421, for example, 
a ^ 0.35. 

Blazar jets are almost certainly in relativistic motion, but although rel- 
ativistic shocks offer a range of indices, the required s — 3.7 is not easily 
obtained. As shown in Fig. 1, mildly relativistic shocks (0.1 < ui < 0.9, 
with ui the speed of the shock seen from the rest frame of the upstream 
plasma) give this value provided the electron pressure is important behind 
the shock front. This effect is due not to a radical change in the way in 
which acceleration operates, but simply to the increased compression ratio 
across the shock front when electrons are heated to a temperature compa- 
rable with their rest-mass. Whether the processes mediating a relativistic 
shock lead to such heating is unknown. However, it seems unlikely that the 
background plasma, which is assumed to carry more energy flux than the 
accelerated component, should remain invisible, despite having an electron 
temperature of several MeV. A more attractive alternative is presented by 
oblique relativistic shocks. In contrast to the nonrelativistic case, where the 
obliquity has no effect on the spectral index (Axford 1980), oblique rela- 
tivistic shocks show much harder spectra (Kirk & Heavens 1989), ranging 
up to s = 3. The reason for this behaviour lies in the increased importance 
of reflections from the shock front as a fundamental accelerating process. If 
particles are tightly bound to held lines, and are also randomly distributed 
in gyro-phase - two assumptions behind the analytic treatment of Kirk & 
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Fig. 1. The power law index of particles accelerated at a relativistic shock front. 
The left-hand panel shows the index s of the phase space-density for a parallel 

shock as a function of the velocity of the shock seen from the rest frame of the 
upstream plasma. Four different equations of state are used: the dashed-dotted 
line is for a relativistic gas, the dotted line for a gas in which the ions provide 
the pressure (i.e., hot ions, cold electrons), the solid line is for hot electrons and 
cold ions, and the dashed line is for a gas containing a population of electron 
positron pairs. A relatively hard spectrum is obtained for intermediate speeds if 
the electrons are hot, or pairs present. At high speeds (F > 5) all curves tend 
towards the value s = 4.2. In the right-hand panel, the effect of a finite angle 
'Pup between the shock normal and the magnetic field is shown, for three different 
values of the upstream shock speed Ush = 0.1, 0.3 and 0.5 and a compression 
ratio of 4. As the speed of the intersection point of the shock and field line (= 
Ush/ cos(#up)) approaches c, the spectra harden. 



Heavens (1989) ~ the probability that they are reflected by the magnetic 
compression at the shock front is significant. As the shock speed increases, 
the energy gained on reflection goes up, since the effective speed of the 'mir- 
ror' is that of the intersection point of a magnetic field line with the shock 
surface. The overall effect is to harden the spectra as shown in Fig. 1. 

It has been pointed out that this type acceleration might be supressed 
in highly relativistic fiows (Bcgelman & Kirk 1990). In the absence of cross- 
field transport, crossing and recrossing requires that a particle move faster 
along a field line than the intersection point of that line with the shock front. 
This speed will normally exceed c for a highly relativistic shock, unless the 
upstream field is aligned to within an angle of l/F of the shock normal. 
(Here F = {1 — tt^/c^)~^/^ is the Lorentz factor associated with the speed 
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of the shock seen from the rest frame of the upstream plasma.) As a result, 
recrossing of the shock into the upstream medium would not be permitted 
and acceleration would be limited to a single encounter - a process usually 
referred to as 'shock-drift' acceleration. 

The two main assumptions used in the analytic theory of relativistic 
oblique shocks (conservation of magnetic moment, and absence of cross- 
field transport) are readily relaxed in a Monte-Carlo simulation, and this 
program has been carried out by several groups (Ostrowski 1991, Lieu et 
al. 1994, Naito & Takahara 1995). The simulations show that the degree of 
cross-field transport is critical. The analytic results are reproduced if it is 
small, but the spectra soften to reach the value for a parallel shock once 
the transport properties become isotropic. Another feature emphasised by 
these simulations is that the spectrum produced depends upon the type of 
scattering used in the simulation - a property also of parallel relativistic 
shocks (Kirk & Schneider 1988). 

Scattering is ultimately just one way of describing the perturbed trajec- 
tory of a particle which encounters fluctuations of the magnetic field about 
its average value. If, as would appear from the above work, the obliquity 
of the shock front is crucial, fluctuations might have an important role in 
changing the angle with which a particular held line hits the shock front, a 
process which is not included in simulations which simply superpose stochas- 
tic scattering on motion in a homogeneous fleld. To investigate this effect, 
two groups have performed simulations in which a tangled magnetic fleld 
with a random component is realised explicitly (Ballard & Heavens 1992, 
Ostrowski 1993). Assuming that the fleld is frozen into the plasma, which 
supports a shock front consisting of a simple velocity discontinuity, the par- 
ticle orbits are followed by numerically integrating the equations of motion. 
Prom the considerations mentioned above, we might expect that reflection 
from the shock front plays the dominant role in such a situation. In fact, a 
tangled fleld will always inhibit the escape of a particle reflected off the ad- 
vancing edge of the shock front and cause it to be caught once again. Thus, 
a layer of energised particles may well be built up on the upstream side of 
the shock front. On transmission to the downstream plasma, tangling of the 
fleld will in principle permit some particles to return to the shock front. In a 
highly oblique shock, this process is unimportant, which may also be the case 
at a relativistic shock with tangled flelds. However, this is mere speculation; 
the two Monte-Carlo simulations unfortunately yield conflicting results, for 
reasons which have not yet been resolved. A possible explanation is that 
the techniques by which the stochastic fleld is realised are different in each 
simulation, and it is not clear that the statistical properties of the resulting 
fleld lines are equivalent. 

Stochastic magnetic flelds have recently been investigated in connection 
with cosmic ray transport (Chuvilgin & Ptuskin 1993) and with nonrela- 
tivistic shocks (Duffy et al. 1995, Kirk, Duffy & Gallant 1996), where they 
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have been shown to be one of the few effects capable of modifying the rather 
robust result given by Eq. (1). Simulation techniques which exploit this ap- 
proach (Gieseler et al. 1997) may prove useful in resolving the question for 
relativistic shocks. 

3 A model of the time dependence of acceleration 

In the case of the homogeneous synchro-self-compton model, the following 
simple estimate can be made (cf. Takahashi et al. 1996). Denote by vio the 
highest energy photons emitted as synchrotron radiation in units of 10 kcV, 
and by the variability timescale of this emission in thousands of seconds, 
which we identify with the cooling time by synchrotron radiation. Then 
vio fa 10~^^7^i?5 and ^3 ^ 10^'y~^ B~'^5~^ , where 7 is the Lorentz factor 
of the particle in the rest frame of the emitting plasma, B the magnetic 
field in gauss and 5 the Doppler boosting factor. These relations imply 



5 = 1 X gauss and 7 = 10^ x tl^^iy^^S''^/^. Those objects 



which show TeV emission generally have roughly the same luminosity in the 
synchrotron and the inverse compton parts of the spectrum, which means 
that the energy densities of photons and magnetic field in the source are 
of comparable magnitude. Setting them equal leads to an expression for 
the 'aspect ratio' of the source region r] = d/R, where d = lOOOt^S/c is 
a characteristic thickness measured in the rest frame of the source, and 
R is defined such that nR^ is the area of the source when projected onto 
the plane of the sky. For Mkn 421, which has an apparent luminosity of 
6 X 10^^ergs-^ one finds rj = 7 x 10-^u^Q^hl^^6^/^, and, inserting 6 = 10, 
vio = 1 and ^3 = 1, we find r] ^ 0.03. This indicates that we are dealing 
with an essentially two-dimensional source region, which it is tempting to 
identify with a relativistic shock front. 

The observation of electrons during their acceleration by the shock front 
of SN1987A (Staveley-Smith et al. 1992) led to the development of a simple 
but useful model for the time-dependence of shock acceleration (Ball & 
Kirk 1992) which can be applied to the present case simply by generalising 
the kinematics to relativistic fiows, and including the effects of synchrotron 
energy losses. The basic idea is to divide the system up into two regions: 
one - which we can imagine to be around or just in front of the shock front 
- in which particles are either repeatedly reflected, or cross and recross the 
front and so undergo continuous acceleration at a rate t"^.^, and a second 
which is the downstream plasma, where there is no further acceleration, but 
only cooling. Escape from the acceleration zone into the downstream plasma 
occurs at a rate t"^. The number A'"(7)d7 of particles in the acceleration 
zone with Lorentz factor between 7 and 7 -|- d7 is governed by the equation 
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(e.g., Kirk, Melrose & Priest 1994), where ps = 4:aT/{Smec){B^/87r) with 
(Tt = 6.65 • 10~^^cm^ the Thomson cross section. The first term in brackets 
in Eq. (2) describes acceleration, the second describes the rate of energy 
loss due to synchrotron radiation averaged over pitch-angle in a magnetic 
field B (in gauss). Particles are assumed to be picked up (injected) into the 
acceleration process with Lorentz factor 70 at a rate Q particles per second. 

In a one-dimensional picture, with the x coordinate along the shock 
normal the kinetic equation governing the differential density d?i(x,7,t) of 
particles which have escaped into the downstream plasma and are in the 
range dx, d7 is 

^-^(ft,^„) = ^.-(x-x.(t„ (3, 

where Xs{t) is the position of the shock front at time t, and a coordinate 
system has been used which in which the plasma is at rest. Note that the 
'injection' term in this equation is provided by those particles escaping the 

acceleration zone. 

Equations (2) and (3) are simple to solve. Using the synchrotron Green's 
function the electron distribution is then readily linked to the emitted radi- 
ation, given the time dependence of the rate Q at which particles are picked 
up by the acceleration process (Kirk, Rieger & Mastichiadis 1997). For ex- 
ample, a simple representation of a flare is found by setting Q{t) = Qq for 
t <0 and t > lOtacc and Q{t) = 2Qo for < t < lOtacc 

The corresponding synchrotron emission in the rest frame of the source 
is shown in Fig. 2. The behaviour of the spectrum at frequencies well be- 
low the maximum, where the acceleration time tacc is much shorter than 
the synchrotron cooling time (/?g7)~^ shows the characteristic 'slow-lag' 
observed in several sources (e.g., Pks 2155-3034: Sembay et al. 1993 and 
Mkn 421: Takahashi et al 1996). This can be interpreted simply in terms of 
synchrotron cooling (Tashiro et al. 1995). However, closer to the maximum 
frequency, the finite acceleration time comes into play (At 7 = 7max) we 
have tacc — (/?s7max)~^-) Then it is possible to produce hysteresis curves 
which are followed in the opposite direction, namely anti-clockwise. This 
behaviour is characteristic of the acceleration mechanism. Although the de- 
tails of such loops will depend on factors neglected here such as the energy 
dependence of taco inhomogeneities in the source geometry, and smoothing 
effects produced by the finite light travel time across the face of the source, 
the conclusion remains that anti-clockwise motion is the signature of the ac- 
celeration mechanism. Combined with the fact that rapid variation implies 
laminar source geometry and that substantial doppler boosting is required, 
this may also be the signature of a relativistic shock front. 
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Fig. 2. The predicted behaviour of the synchrotron emission and spectral index 
during a flare. In the left-hand set of three panels, the light curve, and spectral 

index at low frequency {ly = J^max/20) are shown (with the index a calculated from 
the flux ratio at v and v/lO). The loop in the a vs. intensity plane is traversed 
in the clockwise direction. In the right-hand panels, the same plots are shown at 
v = COi^max. Here the loop is followed in the anti-clockwise direction. 
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